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Abstract. The dynamical response of the solar corona to surface and sub-surface perturbations depends on the chromospheric 
stratification, and specifically on how efficiently these layers reflect or transmit incoming Alfven waves. While it would be 
desirable to include the chromospheric layers in the numerical simulations used to study such phenomena, that is most often 
not feasible. We defined and tested a simple approximation allowing the study of coronal phenomena while taking into account 
a parametrised chromospheric reflectivity. We addressed the problems of the transmission of the surface rotation to the corona 
and that of the generation of polar plumes by Alfven waves [1,2]. We found that a high (yet partial) effective chromospheric 
reflectivity is required to properly describe the angular momentum balance in the corona and the way the surface differential 
rotation is transmitted upwards. Alfven wave-driven polar plumes maintain their properties for a wide range of values for the 
reflectivity, but they become bursty (and eventually disrupt) when the limit of total reflection is attained. 
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1. INTRODUCTION 

The dynamics of the outer layers of the Sun strongly 
constraint the global topology of the solar corona and 
are at the origin of many transient events observed. The 
time- varying magnetic field generated by the solar dy- 
namo connects the sub-surface layers of Sun up to the 
corona and heliosphere, and effectively couples the dy- 
namics of those regions. The photospheric movements 
are expected to shear the coronal magnetic arcades at 
their footpoints. This shearing may be transmitted up- 
wards by the MHD wave-modes and so affect the dy- 
namics of the solar wind and of the corona. The range 
of characteristic time and spatial scales for these sur- 
face motions is very broad. Slow (null frequency) mo- 
tions are widely believed to contribute to the quasi-steady 
accumulation of magnetic energy in the corona, as re- 
quired e.g., by the current storage- and-release type erup- 
tion models. Finite frequency movements, on the other 
hand, inject waves into the same magnetic structures and 
provide for the vast wealth of oscillatory phenomena de- 
tected in the corona. The standard approach to the mod- 
elling of the effects of surface movements on the corona 
is that of the line-tied forcing of the magnetic loops by 
their footpoints, which are rigidly attached to the photo- 
sphere. An important consequence of this approximation 
is that the lower coronal boundary acts as perfect reflec- 
tor to MHD waves instead of letting the chromospheric 
layers define the amount of reflection occurring there. 
Total reflection occurs at the foot-points of the loops, 



which act as resonant cavities. Finite-frequency footpoint 
motions induce long-lasting large amplitude oscillations 
[3] and null frequency motions shear the loops by arbi- 
trarily high amounts. Of man [4] pointed out that foot- 
point leakage of coronal loop oscillations could be im- 
portant in some cases. Grappin et al. [5] have shown that 
the line-tying approximation is inappropriate to charac- 
terise the coronal loop's dynamical response to footpoint 
motions even for very high effective chromospheric re- 
flectivities. An alternative approximation which diamet- 
rically opposes the line-tied condition relies on the so- 
called fully transparent boundaries, as used by Pinto 
et al. [1, 2]. In these works, the authors neglected all 
chromospheric reflection and purposely considered "pes- 
simistic" configurations in what concerns the destabil- 
isation of coronal structures and searched for mecha- 
nisms which nonetheless are able to produce coronal 
events such as polar plumes, jets and inflows. Both ap- 
proximations are extreme ones, the physical reality lay- 
ing necessarily in between. The line-tied case is prone 
to severely overestimate the energy input from the pho- 
tosphere into the corona within magnetically connected 
regions, and the fully transparent case underestimates it 
is the absence of the denser chromospheric (and photo- 
spheric) layers. Ideally, one would include these layers in 
the numerical domain. But that demands prohibitive in- 
creases in CPU time and model complexity, specially in 
multi-dimensional configurations. Here, we test a semi- 
reflective lower boundary condition applied at the base 
of the corona which represents the effective response of 



the chromosphere to impinging Alfven waves. This con- 
dition is comparable to that used by Verdini et al. [6] 
but with a different implementation, more appropriated 
to our MHD setup. 



2. NUMERICAL MODEL 

We used the DIP numerical code, which is a 2.5 D ax- 
isymmetric model of the solar corona solving the fol- 
lowing system of MHD equations describing a one-fluid, 
isothermal, fully ionised and compressible plasma: 
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The magnetic field B decomposes into a time- 
independent external component B° and an induced 
one b. The solar wind develops into a stable transsonic 
solution where the magnetic tension cannot counteract 
its dynamical pressure, hence forming the open-field 
regions of the corona. The upper boundary (placed at 
r = 15) is transparent. The lower boundary (placed 
at r = 1.01 Rq ) is semi-reflective with respect to the 
Alfven mode, but transparent with respect to all others 
(cf. Sect. 2.1). We recur to the description of the system 
of equations above in terms of the MHD characteristics 
at the domain's boundaries in order to achieve total trans- 
parency or to control explicitly its reflectivity. Alfven 
waves are injected at the lower boundary by perturbing 
the corresponding characteristic there. The diffusive 
terms are adapted so that grid scale (Al) fluctuations are 
correctly damped. The kinematic viscosity is defined as 
V = v (A//A/ ) 2 , typically with v = 2 x 10 14 cm 2 • s" 1 
and 0.01 < (A//A/ ) 2 < 10. The grid is non-uniform and 
Al is minimal close to the lower boundary. The magnetic 
diffusivity 77 is scaled similarly. Both coefficients were 
kept as low as possible in all the runs performed here 
(yet inevitably larger than those in the real Sun). Grappin 
et al. [7] give a more thorough description of the spatial 
and temporal schemes used. 



2.1. Semi-reflective lower boundary 

As stated in Sect. 2, our lower and upper boundary 
conditions consist of imposing the amplitudes of the 
MHD characteristics propagating into the numerical do- 
main (the outgoing ones being already completely de- 
termined by the dynamics of the system). We treat the 
chromosphere as a discontinuity, hence the effective re- 
flectivity (of the Alfven mode) is completely determined 



by the ratio £ of Alfven wave speeds above and below. 
The chromospheric reflectivity a is then defined as 
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Ca representing the Alfven speed B/ (/lop) 1 ^ 2 . The ap- 
proximation is valid for perturbations whose characteris- 
tic wavelength is much larger than the thickness of the 
chromosphere. We defined two types of semi-reflective 
conditions, one for the Alfven characteristic and another 
one for purely azimuthal motions (note that, on one hand, 
perturbations with purely azimuthal velocity are a mix- 
ture of MHD wave modes, and on the other hand that 
pure Alfven modes are circularly polarised). The first 
kind translates simply in 
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where, L\ is the amplitude of the upward-propagating 
Alfven characteristic (i.e, propagating into the domain) 
and is the amplitude of its downward-propagating 
counterpart. The second kind assumes the following 
three conditions 
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where f(6,t) describes a perturbation exerted at the sur- 
face. This results in a combination of perturbations of the 
slow-mode, the fast-mode and Alfven mode characteris- 
tics, all defined as a function of f(0,t). The two kinds 
of semi-reflective boundaries are used here with differ- 
ent purposes. We use the first kind in Sect. 4, which ad- 
dresses a problem where the injection and reflection of 
pure and finite-frequency Alfven waves is of importance. 
We favour the second kind of semi-reflective boundary in 
Sect. 3, on which we study quasi-steady solar wind solu- 
tions for which a better control of the surface azimuthal 
velocity is required (and the effects of perturbing non- 
Alfvenic modes are not important). 



3. CORONAL ROTATION 

Pinto et al. [2] studied the response of the solar wind and 
corona to the slow (quasi- steady) variations of the global 
magnetic field during a solar cycle. They produced all 
yr time-series of magnetic fields generated by a 2.5D ax- 
isymmetric kinematic dynamo code [STELEM, 8] which 
was used to constraint the magnetic coronal topology 
in the coronal MHD numerical code DIP. They quanti- 
fied the temporal evolution of the poloidal distribution of 
the wind velocity, mass and angular momentum fluxes, 
which were found to depend strongly on the magnetic 
topology (hence on the moment of the solar cycle). The 



FIGURE 1. Coronal magnetic field, wind velocity and rotation rate at two illustrative moments of the solar cycle (t = and 
t = 3.8 yr in Fig. 3 of [2]). The blue lines are magnetic field-lines. The blue-yellow colour-scale represents the wind speed in Mach 
number (blue: v r /c s < 1; yellow: v r /c s > 1). The orange colour-scale represents the rotation rate £1 The plots on the right side 
represent the rotation period as a function of latitude at different radii (r = 1, 2.5, 7, 15 Rq). 



determination of the angular momentum flux relied on 
the assumption of a solid-body and time-invariant so- 
lar rotation. We revisit here this work and apply semi- 
reflective conditions to the lower boundary of the numer- 
ical domain. We selected two illustrative instants of the 
Pinto et al. [2] simulations. These correspond to t = 
(activity minimum) and t = 3.8 yr (activity maximum), 
and are represented in the first and fourth panels of their 
Fig. 3. We started off from a non-rotating state and ap- 
plied a torque at the surface which accelerated it to the 
following rotation rate profile 

a (6) =Q a + £l b sin 2 6 + a c sin 4 0, (4) 

with 6 being the latitude, Q. a = 14.713 7 da v> &b = 
-2.396 °/day and £l c = -1.787 °/day. Equation (4) fits 
the time-averaged differential rotation pattern at the sur- 
face of the Sun [9]. The duration of the initial accelera- 
tion was ~ 1 /4 of the average (final) rotation period at 
the surface. The initial transient propagates upwards as 
an Alfven wavefront, accelerating the open field plasma 
in the azimuthal direction and exciting a few global os- 
cillations in the closed field regions. After a few Alfven 
transit times, and for a small enough £, the corona and 
wind settle down into a quasi-steady state. For £ ~ 1, the 
highly transparent surface spins down quickly. Note that 
the rotating wind flow carries a net outward angular mo- 
mentum flux, and that a magnetic braking torque is ap- 
plied at the surface in result. For £ ~ 0, the surface main- 
tains the imposed rotation, smaller closed loops keep os- 
cillating resonantly for a long time and the larger ones are 
sheared indefinitely. Intermediate solar-like values allow 
for the surface rotation to be maintained in the open-flux 
regions while the minimal required amount of footpoint 
leakage is allowed for the closed-flux regions to stabilise. 
The streamers gradually assume a nearly rigid rotation, 
while the open field regions (coronal holes) shape up as 
a Parker's spiral. The rotation rate of the former is deter- 
mined by that of its magnetic footpoints. Large stream- 
ers like that on the first panel of Fig. 1 encompass a wide 



latitudinal range, and therefore a wide range of footpoint 
rotation rates. As a result, the surface differential rota- 
tion pattern reflects more visibly in such streamers than 
on the smaller ones. On the overall, the strongest latitudi- 
nal shears appear at the streamer/coronal hole frontiers. 
The finite magnetic resistivity affects the width of these 
regions, but has a negligible effect on the overall rota- 
tion rates. Figure 1 represents the coronal magnetic field 
geometry, the wind speed (in units of Mach number) and 
the coronal rotation rate for t = and t = 3.8 yr. The plots 
to the right of the figures shows the rotation period as a 
function of latitude at different radii (r = 1, 2.5, 7 and 
15 Rq). We note that in order to achieve long-lasting (i.e, 
stable) coronal rotation profiles we had to resort to partic- 
ularly low values for £. The calculations presented here 
were made with £ = 10 -3 rather than a more solar- like 
£ — 10 -2 . The reason for this is that after the initial and 
short-lived forcing, the surface was left to evolve under 
the influence of the braking torque exerted by the wind 
flow alone without any counter-balancing torque due to 
the underlying convective dynamics. This issue will be 
addressed in a forthcoming paper. 

4. POLAR PLUME GENERATION 

Pinto et al. [1] studied the properties of polar plumes gen- 
erated by Alfvenic torsional motions of the footpoints of 
a bipolar magnetic structure in a coronal hole using a 
fully transparent lower boundary. Such magnetic struc- 
tures appear naturally in unipolar regions such as the po- 
lar coronal holes. They are also visible in simulations 
such as that in Fig. 1. We extended this work by account- 
ing for the wave reflection occurring at the dense chro- 
mospheric layers by means of the reflectivity parameter 
£. As before, Alfven waves are continuously injected at 
the coronal base and drive the formation of a dense jet 
(plume) along the magnetic axis of the structure. The 




FIGURE 2. Left panel: the standard case, with £ = 1. Right 
panel: £ = 0.01. The colour-scale represents the density. White 
lines are magnetic field-lines. 
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FIGURE 3. Density n as a function of time at the bottom of 
the jet (r = 0.2 ~ 1 .4 R ). Each curve corresponds to a given £. 
Note the regime change happening for £ = 0.01. The ordinates 
are in log-scale. 



density and velocity contrast between the axis of the 
plume and the background wind flow decays with height, 
unlike in thermally driven plume models [10]. The jet 
displays a series of "blobs" (not to be confused with 
streamer blobs [11]) propagating outwards along its axis, 
which are mostly visible in the first 7 solar radii and can 
be identified as slow mode wave-fronts generated non- 
linearly at the base of the plume [cf. 12, 13]. The growth 
rate now depends both on the frequency of the surface 
motions and the effective chromospheric reflectivity. Fig- 
ure 2 shows a snapshot of the density at t = 7.5 h after the 
beginning of the wave injection for the case with £ = 1 
(the standard case with transparent lower boundary) and 
£ = 10 -2 . The colour-scale represents the density cor- 
rected by a factor (r/ R Q ) 2 x exp (1/ R Q — 1/r) to com- 
pensate for the background stratification and let the den- 
sity structures be visible. The white-lines are magnetic 
field- lines. Figure 3 shows the temporal evolution of the 
density at the axis of the plume, just above the null point. 
Each curve corresponds to a different £. The system's 
evolution remains qualitatively similar for 0.1 < £ < 1, 
but changes dramatically for £ < 0.1. The flow along the 



axis switches to a regime where it is dominated by the 
slow-mode wave-fronts (the "blobs") rather than by the 
smooth overdense upflow. For £ w 0, the plumes are eas- 
ily disrupted. 

5. CONCLUSIONS 

We studied two different problems concerning the cou- 
pling between surface movements and the solar corona. 
The first one consists of the transmission of the sur- 
face rotation to the corona and the second one of the 
generation of polar plumes by Alfven waves. For this 
purpose, we implemented and tested a semi-reflective 
lower boundary condition for our coronal model which 
accounts for the Alfven wave partial reflection taking 
place in the chromosphere. The effective chromospheric 
reflectivity is parametrised by the ratio of Alfven speeds 
across the chromosphere and transition region. We found 
that a high (yet partial) effective reflectivity is required 
for sustaining the coronal rotation against the solar wind 
magnetic breaking torque, while still allowing for the 
necessary amount of footpoint leakage. Future work will 
consider the underlying time-dependent convective sur- 
face dynamics. We verified that the generation of polar 
plumes by Alfven waves works well beyond the initially 
proposed limit of full transparency [1]. The plume prop- 
erties are maintained for a wide range of values for the 
reflectivity, but become bursty (and eventually disrupt) 
when the limit of total reflection is reached. Additional 
work is in progress to evaluate more finely the parameter 
range < £ < 10~ 2 . 
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